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Abstract. - The momentum and energy dependence of the weight distribution in the vicinity 
of the one-electron spectral-function singular branch lines of the ID Hubbard model is studied 
for all values of the electronic density and on-site repulsion U. To achieve this goal we use 
the recently introduced pseudofermion dynamical theory. Our predictions agree quantitatively 
for the whole momentum and energy bandwidth with the peak dispersions observed by angle- 
resolved photoelectron spectroscopy in the quasi-lD organic conductor TTF-TCNQ. 



The finite-energy spectral dispersions recently observed in quasi-one-dimensional (ID) met- 
als by angle-resolved photoelectron spectroscopy (ARPES) reveal significant discrepancies 
from the conventional band-structure description [1,2]. The study of the microscopic mecha- 
nisms behind these unusual finite-energy spectral properties remains until now an interesting 
open problem. There is some evidence that the correlation effects described by the ID Hub- 
bard model might contain such finite-energy mechanisms [1,2]. However, for finite values 
of the on-site repulsion U very little is known about its finite-energy spectral properties, in 
contrast to simpler models [3] . Bosonization [4] and conformal-field theory [5] do not apply at 
finite energy. For U — ► oo the method of Ref. [6] provides valuable qualitative information, yet 
a quantitative description of the finite-energy spectral properties of quasi-lD metals requires 
the solution of the problem for finite values of U. The method of Ref. [7] refers to features of 
the insulator phase. For U ~ 4i, where t is the transfer integral, there are numerical results 
for the one-electron spectral function [8] which, unfortunately, provide very little information 
about the microscopic mechanisms behind the finite-energy spectral properties. Recent pre- 
liminary results obtained by use of the finite-energy holon and spinon description introduced 
in Refs. [9-11] predict separate one-electron charge and spin spectral branch lines [1]. For 
the electron-removal spectral function these lines show quantitative agreement with the peak 
dispersions observed by ARPES in the quasi-lD organic conductor TTF-TCNQ [1]. However, 
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these preliminary studies provide no information about the momentum and energy depen- 
dence of the weight distribution in the vicinity of the charge and spin branch lines and do not 
describe the TTF dispersion. The main goal of this paper is the evaluation of such a depen- 
dence for all values of U and electronic density. In order to solve this complex many-electron 
problem, we use the pseudofermion dynamical theory recently introduced in Ref. [12]. 

The model reads H = —t^lj A c ] <j c j+i,° + ^- c -] + %',T^>i where a and Cj t „ 

are spin a =T,| electron operators at site j = l,...,N a and h^ a = cj a Cj, a . The low- 
energy spectral properties of TTF-TCNQ involve inter-chain hopping and electron-phonon 
interactions. Thus, our results are to be applied above the energies of these processes. We 
consider an electronic density n — N/N a in the range < n < 1 and zero magnetization 
where N is the electron number. The Fermi momentum is kp — nn/2 and the electronic 
charge reads — e. The one-electron spectral function B l (k, ui) such that Z = — 1 (and / = +1) 
for electron removal (and addition) reads, B^ 1 (k, u) = J2 a 7 I (tI c k,a\ GS)\ 2 S(uj + AE^ 1 ) 

and B +1 (k, u) = Ha, 1 '\ l a'\ c l^\ GS )\ 25 ^ ~ AE* +1 ). Here c Ka and c\ a are electron 
operators of momentum k and \GS) denotes the initial A-clcctron ground state. The 7 and 7' 
summations run over the N — 1 and N + 1-clectron excited states, respectively, and AE^ 1 
and AE* +1 are the corresponding excitation energies. 

The pseudofermion dynamical theory reveals and characterizes the dominant microscopic 
processes which generate over 99% of the electronic weight of the functions B l (k, lu), and also 
describes all other processes [12]. The weight distribution in the vicinity of the branch lines we 
consider below is fully controlled by such dominant processes. Those involve the —1/2 Yang 
holons and the cO, si, and cl pscudofcrmions studied in Ref. [11]. For simplicity, we denote the 
first two objects by c and s pseudofermions, respectively. The c pseudofermion carries charge 
— e and has no spin and the s pseudofermion is a spin-zero two-spinon composite object and has 
no charge. The cl pseudofermion is a ?7-spin-zero two-holon composite object, carries charge 
— 2e, and has zero spin. The —1/2 Yang holon has r?-spin 1/2, 77-spin projection —1/2, charge 
— 2e, and zero spin. The c, s, and cl pseudofermions carry momentum q = q + Q a (q)/N a , 
where in our case a = c, s, cl. Here q is the bare-momentum and Q a (q)/N a is a momentum 
functional defined in Ref. [11], whose expression involves the two-pseudofermion phase shifts 
& aa '(q,q') defined in the same reference, where a, a' = c,s,cl. Following the one-to-one 
correspondence between the momentum q and bare-momentum q, one can either label a a 
pseudofermion by q or q. Here we use the bare-momentum q. The above pseudofermions 
have energy bands e c (q), e s (q), and e cl {q) = E u + e° cl {q) such that \ q\ < n, \q\ < kp, and 
I q\ < I 71 " — 2fcp], respectively. These bands arc studied and plotted in Ref. [10]. E u defines the 
lower-limit of the upper- Hubbard band (UHB) [9] and equals the energy required for creation 
of a —1/2 Yang holon, which is a dispersion-less object. It is such that E u = 4tcos(7rn/2) 
for U/t — » 0, E u = U + 4£cos(7rn) for U » t, E u = U + At for n — > 0, and as n — > 1 E u 
approaches the value of the Mott-Hubbard gap [10]. In the ground state there are no —1/2 
Yang holons, the cl and s pseudofermions bands are empty and filled, respectively, and the c 
pseudofermions occupy < | q\ < 2kp (leaving 2kp < \ q\ < n empty). 

The ground state and excited states can be expressed in terms of occupancy configurations 
of the above quantum objects. For electron removal, the dominant processes involve creation 
of one hole both in e c (q) and e s (q)- For electron addition, these dominant processes lead to two 
structures: A lower- Hubbard band (LHB) generated by creation of one particle in e c (q) and one 
hole in e s (q); A UHB generated by creation of one hole both in e c (q) and e s (q) and either one 
particle in f° cl (q) for n < 1 or one —1/2 Yang holon for n — > 1. According to the pseudofermion 
dynamical theory of Ref. [12], both the one-electron spectral- weight singularities and edges 
are located on pseudofermion branch lines. Such lines are generated by processes where a 
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Fig. 1 - Extended-zone scheme centered at k = for U — 4.9 1 and n = 0.59. The solid and dashed 
lines denoted by the letters c, c', c", and s are singular and edge branch lines, respectively. Electron 
removal (LHB addition) corresponds to ui < (and ui > 0) and to — E u marks the UHB lower limit. 



specific pseudofcrmion is created or annihilated for the available values of bare-momentum 
q and the remaining quantum objects arc created or annihilated at their Fermi points. The 
weight shape of the singular (and edge) branch lines is controlled by negative (and positive) 
exponents smaller than zero (and one). The electron removal (to < 0) and LHB addition 
(u) > 0) singular and edge branch lines are represented in Fig. 1 by solid and dashed lines, 
respectively. For simplicity, the figure docs not represent the u) > E u UHB region. (Below 
we find that the shape of the UHB singular branch lines is fully determined by the shape of 
the electron-removal singular branch lines represented in Fig. 1.) The dashed-dotted lines 
and some of the branch lines of Fig. 1 are border lines for the u < E u domain of the (k, em- 
plane whose spectral weight is generated by dominant processes. (There is a region limited 
above by the s line for kp < k < 3k F and below by the c" and c' lines for k F < k < 2kp and 
2k F < k < 3k F , respectively, which does not belong to that domain.) The dominant processes 
also include particle-hole pseudofcrmion processes which lead to spectral weight both inside 
and outside but in the close vicinity of that domain. 

For simplicity, in this letter we consider the exponents that control the weight shape in 
the vicinity of singular branch lines only. However, similar expressions can be evaluated for 
any branch line. The singular branch lines correspond to all the (fc, w)-plane regions where 
there are weight-distribution singularities. While it is difficult to measure the exponents 
experimentally, a crucial test for the suitability of the model to describe real quasi- ID materials 
is whether the ARPES peak dispersions correspond to the predicted singular branch lines. The 
general branch-line spectral-function expression given in Rcf. [12] applies provided that the 
excited states associated with the one-electron branch lines are used. 

We start by considering the spin s branch line for < k < 3k F , the charge c branch line 
for < k < 7T — kp, and the charge c' branch line for < k < it + kp (see Fig. 1). The 
parametric equations that define these branch lines read u(k) = e s (q) for the s line where 
q = q(k) = (l + l) kp-lk for (1 + l)k F /2 < k < kp + (1 + l)k F and uj(k) = e a (q) for the a = c 
line (i' = +1) and a = c' line (i' = —1) where q = q(k) = k + l 1 kp for < k < n — l 1 kp. 
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Here, I — ±1 and e c >(g) = e c (q)- The following expression describes the weight distribution 
in the vicinity of the a = s, c, d branch lines for lo values such that (e a (q) + Jw) is small and 
positive, 

B l {k, w) = C l a (k) (e a {q) + Ilo) " ; a = s,c,c'. (1) 

The ^-dependence of C l a (k), such that C l a (k) > for U/t > 0, is in general involved and can 
be studied numcr ically. For U/t 0, C l a (k) behaves as C l s {k) Si _ x C s , C^(fc) -► <5 i;+1 C* c , 
and C l c ,(k) — ► 0, where C s and C c are independent of k. The exponent of Eq. JIJ reads, 



C.(*0 



-i+ y [j_L- A $ sa 



2£o 4 



14, 



where a = c for t/ = +1, a = c' for tl = —1, and fi(x) = 1 — / (1 + x). In equation 
the phase shifts are defined in Ref. [11] and £o = \J^K p where K p , such that K p — > 1 as 
U/t —> and -Kp — > 1/2 as J7/t — » oo, is defined in Ref. [4]. The exponents ( s (k) and 
Ca=c,c' (k) are plotted in Figs. 2 and 3, respectively, as a function of fc for several values of 
U/t and n = 0.59. The exponent ( c {k) (and Cc'(fc)) is plotted in Fig. 3 (a) and (b) (and Fig. 
3 (c) and (d)) for I = — 1 and I = +1, respectively. As i//t — > oo these exponents behave 
as C s (fe) = -1/2 + 2(q/4fc F ) 2 and Cs(fc) = -(g/2fc F )[l - (q/4k F )} for i = -1 and I = +1, 
respectively, and Cc(s) = Cc'(<z) = — 3/8, in agreement with Ref. [6]. The weight shape of the 
UHB singular branch lines is controlled by the same exponents as the corresponding electron- 
removal branch lines. There is a UHB s u branch line for n — k F < k < tt, c u branch line for 
7r — kp < k < 7r, and c' u branch line for n — 3k F < k < tt. For < k < tt and u> > E u the 
parametric equations of these lines is the same as for the corresponding s, c, and d branch 
lines, respectively, provided that k is replaced by it — k and to by E u — u> in the parametric 
equations of the latter lines. Under such a replacement the exponents that control the weight 
shape of these UHB lines are the same as for the corresponding electron-removal lines. While 
in the limit n — > 1 such a correspondence refers also to the value of the constant C~ 1 (fc) of 
Eq. (JJJ, otherwise the corresponding UHB constant C^ 1 (7r — k) is slightly smaller. This is 
because for n < 1 there is also weight in the vicinity of four UHB cl pseudofermion branch 
lines. We omit here the study of the weight shape of these lines whose exponents are positive. 

When the exponents (J2J tend to — 1 as U/t — > the expression is replaced by S(e a (q) + 
I lo). For < k < tt and U/t — > all spectral weight is transferred over to the s branch line for 
< k < kp, c branch line for kp < k < tt — kp, and s u branch line for tt — kp < k < tt. There 
is no weight in other branch lines and k > regions of the (k, a->)-plane in such a limit. The 
two points (k = tt — kp, uj = E u ) and (k = tt — kp, lo = e c (Tr)) where q — tt corresponds to 
k = TT—kp become the same point as U/t — > and the c and s u branch lines become connected 
at that point. By use of the U/t — * expressions given in Ref. [10] for e c (q) and e s (q), one finds 
that these three branch lines give rise to the electronic spectrum u>(k) = — 2i[cos(fc) — cos(k F )]. 
Consistently, the above corresponding exponents are such that Cs(&) — * —1 for < fc < kp, 
Cc(k) — > —1 for kp < k < tt — kp, and Cs u (k) = Cs( 7r — k) — > —1 for tt — kp < k < tt as 
U/t — > 0. Then the correct non-interacting electronic spectral function is reached as U/t — > 0. 
(For U/t — > the exponents 10 behave as Cs(^) = I — ±U Ca{q) = —l/V% + 1/2 for a = c, d 
and I = —1, and ( c {q) = —1 and ( c '(q) = 1 for / = +1.) 
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Fig. 2 - Momentum dependence along the spin s branch line of the exponent ©• 

We note that for (k, u>) points in the vicinity of the branch-line end points (fco, jE u ) 
where J = for a = s,c,c' and j = 1 for a = s u ,c u ,c' u , respectively, the spectral-function 
expression remains valid provided that the exponent £ a [k) given in Eq. J5J is replaced 
by Q(k) = Ca(k) - 2Af(fc) for uj » -le a (q(k)) = ±lv a (k - k ) + jE u . Here I = ±1, 
v a = de a (q)/dq\ q=qFa , qp c — 2fc_F, q Fs — kp, and 2AJ 1 (fc) is a functional denned in Ref. [12]. 
(In Fig. 1, k = kp, k — 3kp, and fc = 5kp for j = 0.) For j = this limit corresponds to the 
so called low-energy Tomanaga-Luttinger liquid regime and the above exponent Ca(^) equals 
that given in Eq. (5.7) of Ref. [5]. The branch-line k domain corresponding to that exponent 
has zero measure relative to the k domain of the branch- line exponent Ca{k)- Therefore, in 
this paper we do not consider such a limiting regime. Another case of the interest is the 
behavior of the spectral function in the vicinity of the end points (fco, jE u ) when these points 
are approached by lines which do not cross the a branch lines. In this case the weight- 
distribution is controlled by an exponent different from both C a (fc) and [13]. 

An interesting realization of a ID metal is the organic charge-transfer salt TTF-TCNQ [2]. 
The experimental dispersions in the electron removal spectrum of this quasi-lD conductor as 
measured by ARPES are shown in Fig. 4. The data were taken with He I radiation (21.2 
eV) at a sample temperature of 60 K on a clean surface obtained by in situ cleavage of a 
single crystal. Instrumental energy and momentum resolution amount to 70 meV and 0.07 
A -1 , respectively. The experimental data in Fig. 4 reproduce earlier work [1,2]. While our 
above theoretical weight-distribution expressions refer to all values of U jt and n, we find 
that the electron removal spectra calculated for t = 0.4 eV, U — 1.96 eV (U/t — 0.49), and 
n = 0.59 yields an almost perfect agreement with the three TCNQ experimental dispersions. 
The exception is the low-energy behavior, as a result of the inter-chain hopping and electron- 
phonon interactions, as mentioned above. If accounted for a renormalization of the transfer 
integral due to a possible surface relaxation [1], these values are in good agreement with 
estimates from other experiments [2] . The experimental TCNQ finite-energy peak dispersions 
of Fig. 4 correspond to the spin s branch line and charge c and c 1 branch lines. Importantly, 
only these singular branch lines, whose weight shape is controlled by negative exponents, lead 
to peak dispersions in the real experiment. The other peak dispersion of Fig. 4 is associated 
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with the electron-removal spectrum of the TTF chains whose density is n — 2 — 0.59 = 1.41. 
It can be described by the electron-addition spectrum of a corresponding n — 0.59 problem. 
Remarkably, we find again that the theoretical lines match the TTF dispersion provided that 
t w 0.27 eV and U < 0.2t within experimental uncertainty. Indeed, for such small values 
oiU/t and n = 0.59 both the electron-addition exponents of Figs. 2 and 3 (d) are positive, 
whereas the exponent of Fig. 3 (b) and the UHB exponent Cs u (k) = Cs(^ — k) for tt — k < hp 
(see Fig. 2) are negative. For n = 1.41 the c and s u lines associated with these exponents 
correspond to electron removal and thus appear at negative values of u> and lead to a single 
singular branch line formed by the c line for 0.597r/2 < k < lAlw/2 and the s u line for 
1.4l7r/2 < k < tt. The energy pseudogap at k — 1.4l7r/2, which separates the c line from 
the s u line, nearly vanishes for U/t < 0.2t and thus is not observed in the experiment. This 
singular branch line is that matching the TTF peak dispersion of Fig. 4. (It matches such 
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Fig. 4 - Angle-resolved photo-emission spectra of TTF-TCNQ measured along the easy transport 
axis and matching theoretical branch lines. 

a dispersion for 0.597r/2 < k < ir.) We thus conclude that in contrast to the conventional 
band structure description, the singular branch lines obtained from the ID Hubbard model 
describe quantitatively, for the whole finite-energy band width, the peak dispersions observed 
by ARPES in TTF-TCNQ. This seems to indicate that the dominant non-perturbative many- 
electron microscopic processes described by the pseudofermion dynamical theory of Ref. [12], 
which control the weight distribution of these charge and spin singular branch lines, also 
control the unusual finite-energy one-electron spectral properties of TTF-TCNQ. 

We thank E. Jeckelmann for stimulating discussions. K.P. thanks the financial support 
of OTKA grants D032689 and T037451, L.M.M. of FCT grant BD/3797/94, and R.C., M.S., 
and U.S. of Deutsche Forschungsgemeinschaft (CL 124/3-3 and SFB 484). 

REFERENCES 

[1] Sing M., Schwingenschlogl U., Claessen R., Blaha P., Carmelo J. M. P., Martelo L. 

M., Sacramento P. D., Dressel M. and Jacobsen C. S., Phys. Rev. B, 68 (2003) 125111. 
[2] Claessen R., Sing M., Schwingenschlogl U., Blaha P., Dressel M. and Jacobsen C. 

S., Phys. Rev. Lett, 88 (2002) 096402; Zwick F., Jerome D., Margaritondo C, Onellion 



8 



EUROPHYSICS LETTERS 



M., Voit J. and Grioni M.,Phys. Rev. Lett, 81 (1998) 2974; Kagoshima S., Nagasawa H. 

and Sambongi T., One-dimensional conductors (Springer, Berlin)1987 and references therein. 
[3] Arikawa M., Saiga Y. and Kuramoto Y., Phys. Rev. Lett, 86 (2001) 3096; Penc K. and 

Shastry B. S., Phys. Rev. B, 65 (2002) 155110. 
[4] Schulz H. J., Phys. Rev. Lett, 64 (1990) 2831. 

[5] Frahm H. and Korepin V. E., Phys. Rev. B, 43 (1991) 5653 and references therein. 

[6] Peng K., Hallberg K., Mila F. and Shiba H., Phys. Rev. Lett, 77 (1996) 1390; ibid. Phys. 

Rev. B, 55 (1997) 15475. 
[7] Sorella S. and Parola A., Phys. Rev. Lett, 76 (1996) 4604. 

[8] Senechal D., Perez D. and Pioro-Ladriere M., Phys. Rev. Lett., 84 (2000) 522. 

[9] Carmelo J. M. P., Roman J. M. and Penc K., Nucl. Phys. B, 683 (2004) 387. 
[10] Carmelo J. M. P. and Sacramento P . D., Phys. Rev. B, 68 (2003) 085104. 

[11] Carmelo J. M. P . , |cond-mat /0305568| 

[12] Carmelo J. M. P. and Penc K., |cond-mat/0311075| 

[13] Carmelo J. M. P., Martelo L. M. and Sacramento P. D., J. Physics: Cond. Matt, 16 
(2004) 1375. 



